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Classically  the  variational  equations  for  the  problem  cf  two  bodies 
have  been  studied  in  Jacobi's  elliptic  elemerts  or  Delaunay's  canonical 
variables  which  make  them  trivial  to  solve.  This  drastic  simplification 
is  achieved  by  successive  canonical  transformations.  Only  recently  astron¬ 
omers  considered  using  rectangular  coordinates  in  Perturbation  Theory, 
which  led  them  to  compute  the  partial  derivatives  of  the  coordinates  along 
a  Keplerian  motion  with  respect  to  Jacobi's  elliptic  elements  (Eckert  and 
Brouwer  1937,  Brouwer  and  Clc.mence  1961).  From  the  Eckert-Brouwer  formula 
for  differential  corrections  cf  Cartesian  coordinates,  Danby  (1965)  de¬ 
rived  a  decomposition  of  the  matrizant  R(t;ty)  into  the  matrix  product 
R(t;tg)  =  M(t)  o  M  *(tg):  the  matrix  M(t)  is  surprisingly  simple  in  ap¬ 
pearance.  but  no  explicit  analytical  expression  is  given  for  its  inverse. 

Similar  to  the  Eckert-Bronwer  formula  is  Bower's  approach  to  the  con¬ 
struction  of  the  matrizant,  (Bower  1932).  It  proceeds  in  two  steps:  the 
variations  of  the  elliptical  elements  with  respect  to  the  initial  conditions 
are  first  computed,  then  the  variations  of  the  coordinates  with  respect  to 
the  elliptical  elements.  The  algorithm  is  left  in  the  form  of  a  sequence 
of  operations  to  be  carried  on  a  desk  calculator  to  arrive  at  the  value  of 
the  matrizant  at  a  given  instant.  Following  Bower's  method,  Danby  (19b2,1965) 
derived  explicit  expressions  for  the  elements  of  the  matrizant  in  the  apsidal 
frame  of  reference  (the  axis  Oz  is  normal  to  the  orbital  plane,  while  Ox 
points  toward  the  periastron  ) .  Such  formulas  had  already  been  produced  by 
Myachin  (1959)  apparently  without  much  explanation  as  to  how  he  obtained  them. 
The  matrizant  in  the  apsidal  coordinate  system  looks  slightly  simpler  when 
the  initial  time  ty  is  taken  at  a  passage  at  periastron  (Danby  1964). 
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Starting  from  the  series  f  and  g,  Sconzo  (1963)  arrives  at  a 
completely  explicit  expression  of  the  matrizant  but,  in  view  of  the  forms 
we  have  given  to  it,  it  is  felt  that  Sconzo's  formulas  could  be  presented 
in  a  more  direct  and  simpler  form. 

Brumberg  (1961)  chooses  to  operate  in  the  orbital  coordinate  system 
(the  axis  Ox  points  towards  the  instantaneous  position).  He  is  as  la¬ 
conic  as  Myachin  about  the  way  he  derives  a  fundamental  set  of  solutions 
to  the  variational  equations  in  this  moving  frame  of  reference.  Charnyi 
(1963)  made  a  first  attempt  at  clarifying  this  issue;  lately  he  has  come 
with  a  rule  to  construct  a  fundamental  set  of  solutions  to  the  adjoint 
variational  equations  for  a  conservative  Hamiltonian  system  with  n  degrees 
of  freedom  when  (n-1)  independent  integrals  are  known  (Charnyi  1965). 

The  author  is  apparently  unaware  of  the  relation  between  his  method  and 
Jacobi's  dual  theorem  of  the  last  multiplier. 

To  conclude,  in  treating  the  variational  equations  of  the  problem  of 
two  bodies,  it  is  felt  that  a  constructive  and  systematic  approach  is  needed 
which  would  permit  an  easy  transformation  from  one  frame  of  reference  to 
another  and  would  also  facilitate  a  proper  choice  of  elements  under  different 
circumstances. 

In  order  to  have  a  more  organized  view  and  to  facilitate  the  search 
for  new  forms  of  the  matrizant,  an  application  of  Jacobi's  dual  theotem  of 
the  last  multiplier  is  made  to  the  solution  of  the  adjoint  variational  equa¬ 
tions  in  an  arbitrarily  given  inertial  coordinate  system.  This  method  per¬ 
mits  the  formation  of  the  matrizant  in  a  straightforward  manner,  and  its 
transformation  to  any  other  coordinate  system  such  as  the  orbital  frame  con¬ 
sidered  by  Brumberg  or  the  intrinsic  frame  postulated  by  Hill's  equation. 
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Only  the  two-dimensional  problem  of  two  bodies  is  examined  in  this 
paper.  Moreover  we  make  quite  apparent  the  places  where  the  treatment  we 
present  exclude  successively  the  rectilinear  solutions,  then  the  circular 
motions  and  finally  the  parabolic  orbits.  These  restrictions  do  not  come 
from  the  introduction  of  a  new  independent  variable  such  as  the  eccentric 
anomaly  and  its  relation  to  the  time  by  means  of  Kepler’s  equation;  rather 
they  follow  from  the  conditions  under  which  the  three  variational  integrals 
are  independent  -,'nd  the  fourth  solution  associated  with  the  last  multiplier 
is  independent  fro;.'  the  gradients  of  the  variational  integrals. 

A  special  effort  is  made  to  keep  the  matrizants  in  the  form  most  suit¬ 
able  to  repeated  processing  by  electronic  computers.  As  far  as  possible  we 
eschew  introducing  transcendental  functions;  the  time  is  kept  as  the  running 
variable  along  the  orbit  so  that  there  is  no  need  for  solving  Kepler’s  equa¬ 
tion.  We  hope  that  formulas  of  this  kind  would  prove  to  be  elementary  enough 
to  fit  into  an  airborne  computer,  when  they  are  extended  to  the  three  dimen¬ 
sional  problem. 

1 •  Jacobi’s  theorem  of  the  last  multiplier  applied  to  the  variational  equa¬ 
tions  . 

As  it  is  usually  done,  we  use  the  conservation  of  linear  momentum  to 
reduce  the  problem  of  two  bodies  to  a  dynamical  system  described  by  the 
Hamiltonian  function 

=  ~  (X2  +  Y2)  -  ^  (1) 

in  a  four  dimensional  phase  space  product  of  the  Cartesian  plane  of  con¬ 
figurations  (x,y)  by  the  Cartesian  plane  of  moments  (X,Y).  In  the  func¬ 
tion  (1),  u  is  a  strictly  positive  real  constant  and 
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r  =  |  x2  +  y2|?2  (2) 

The  Hamiltonian  equations  generated  from  (1)  will  be  called  orbital.  They 
admit  four  integrals: 

(i)  the  energy 

H  =  y  (X2  +  Y2)  -  ^  ;  (3) 

(ii)  the  angular  momentum 

G  =  xY  -  yX;  (4) 

(iii)  Laplaca's  vector  given  by  its  Cartesisn  components 

P  =  GY  -  u~  (3) 

Q  =  -  CX  -  iP;  (6) 

The  functions  (3)-(6)  will  be  termed  orbital  integrals .  We  agree  to  denote 
by  the  same  symbol  either  the  function  which  is  an  integral  or  the  constant 
value  which  this  function  takes  along  a  particular  solution  of  the  orbital 
equations;  no  confusion  could  result  from  this  ambiguity  in  the  notations. 

We  recall  an  essential  formula  relating  the  four  orbital  integrals, 
namely 

P2  +  Q2  =  u2  +  2HG2,  (7) 

and  we  set  up  a  list  of  basic  identities 

2 

Px  +  Qy  =  G  -  ur, 

Py  _  Ox  =  Grr, 

PX  +  QY  =  -ur, 

PY  -  OX  =  G(V2  -  ^) 


(8) 


which  we  shall  use  quite  often  in  the  following  sections  of  this  paper 
without  explicitly  referring  the  reader  to  them. 


Let 


f:  t  ■*  (x(t),  y(t),  X(t),  Y(  t) ) 


(9) 


be  a  solution  of 
variation  of  I' 
tional  equations 
u  =  U 

w  =  W 


the  orbital  equations.  By  definition  a  displacement  or 
is  any  solution  u,  =  (u,w,U,W)  of  the  so-called  varia- 


(10) 


considered  along  the  orbit  P. 


The  variational  equations  are  Hamiltonian:  they  are  generated  from  the 
Hamiltonian  function 


12  2 
0/=  j(tr  +  VO 


(ID 


which,  as  it  is  expected,  is  a  quadratic  form  in  the  vector  u,  whose  coef¬ 
ficients  are  functions  of  the  time  to  be  computed  along  the  orbit  (9). 

Each  orbital  integral  F  generates  an  integral  <5F  of  the  variational 
equations,  which  we  shall  call  here  a  variational  integral.  It  is  to  be  built 
as  the  scalar  product  <$F  =  (uJvF)  of  the  vector  £  by  the  gradient 


VF  -  (F  ,  F  ,  Fy,  Fy)  of  the  function  F  in  the  phase  space  (x,  y,  X,  Y). 

x  y  * 
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Let  us  compute  the  gradients  VC,  VP  and  VQ  of  the  orbital  inte¬ 
grals  (4),  (5)  and  (6): 


G=Y.  p  =  -  u 

x  *  x  3 

r 


G  =  -X,  P  =  -XY  +  y 

y  y  3 

r 


gx  =  -y,  px  -  -W 


Qx  -  -XY  +  u  & 
r 

n  -  y2  x2 
Q  -  X  -  u  — r 

y  r3 


Qx  =  2yX  -  XY 


Gy  -  X,  Py  =  2xY  -  yX 


qy  =  -XX. 


We  shall  now  build  the  six-dimensional  2-vector  VP  a  VQ  which  is  the  ex¬ 
terior  product  of  the  vector 


by  the  vector 


tp  ■  Px  %  +  Py  %  +  Px  %  +  PY  *4 


?Q  *  Qx  %  +  Qy  %  +  Qx  %  +  QY  V 


The  components  of  the  result  written  as  the  sum 

VP  A  VQ  =  2  Aij  %  A  ^ 

(i.j)  J  J 

are  found  to  be 


12  = 

3 

r 

2X 

A23 

13  = 

-  G  |y2  - 

CM 

< 

>— * 

(1 

G^XY  - 

2u^) 

A34 

;(XY  -  2“  f) 


in  the  same  manner  we  construct  the  exterior  products 


VG  A  VP 


.  zL  &i-  A  £-i 


<i,J) 


VG  A  VQ  =  7  C..e.Ae. 

(1.3) 


whose  components  are  as  follows: 


B12 

3 

-MC-i 

C12 

= 

-UG 

r 

r 

3 

2 

B13 

= 

-u  *3 

C13 

= 

-GY  +  u 

3 

r 

2 

O 

i. 

r 

bia 

= 

GY  +  v 

SL. 

3 

C14 

= 

x  y 

-u  -f 

r 

r 

2 

2 

B23 

= 

xy 

U  3 

C23 

= 

GX  -  p 

2LX 

3 

r 

2 

3 

r 

B2A 

s 

-GX  -  1 

.  2LJL 

J  3 

C2A 

= 

X 

u  — 

r 

r 

B3A 

= 

-Gy 

C3A 

= 

Gx. 

(14) 


The  following  two  properties  are  equivalent: 

a)  The  gradients  VP,  VQ  and  VG  are  linearly  independent ; 

£>)  The  Keplerian  motion  r  is  not  rectilinear. 

Indeed  let  us  build  the  four-dimensional  3-vector 

VG  A  VP  a  VQ  -  \  D,  e 


(i 


.j.k)  ijk  **  A  **  A 


Its  coefficients  are  as  follows 


°123  = 

-UG2 

JL 

3 

D134 

r 

D124  “ 

UG2 

JL. 

3 

U234 

=  -g2x. 


(15) 
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But  the  gradients  VP,  VQ  and  VG  are  linearly  dependent  if  and  only  if 
the  3-vector  VG  a  A  VQ  vanishes  identically.  As  it  is  exhibited  by 


the  components  found  in  (15),  this  is  the  case  if  and  only  if  the 

angular  momentum  G  is  zero,  thus  if  and  only  if  the  Keplerian  motion  is 


rectilinear. 


From  now  on  in  this  paper,  we  shall  assume  that  G  is  a  nonvanishing 
constant;  provided  the  plane  (x,y)  of  configurations  is  given  the  good 
orientation,  we  car.  even  assume  that  G  is  strictly  positive. 

Since  G  is  ^  0,  the  fundamental  identity  (7)  implies  that 

7H  «  -  VG  +  ~  VP  +  VQ  (16) 

G  G2  G2 

which  means  that  the  gradient  VH  of  the  orbital  energy  belongs  to  the  three 
dimensional  vector  subspace  generated  by  the  gradients  VG,  VP  and  VQ.  Con 
sequently,  the  orbital  integrals  produce  only  three  linearly  independent  in¬ 
tegrals,  namely  the  functions: 

6G  =  (ii|  VC)  =  Gxu  +  Gyw  +  GXU  +  G^W, 

<SP  =  (u|7P)  =  P^u  +  Pyw  +  PXU  +  PYW,  (17) 

<$Q  =  (uJVQ)  =  Qxu  +  QyW  +  QXU  +  QyW. 

Now,  because  they  are  Hamiltonian,  the  variational  equations  have  a 

multiplier,  which  is  equal  to  1.  Accordingly,  in  view  of  Jacobi’s  theorem 
of  the  last  multiplier,  the  variational  equations  are  to  be  solved  by  quad¬ 


ratures  only. 
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Like  the  variational  equations,  the  adjoint  variational  equations  are  Hamil¬ 
tonian;  indeed  they  are  generated  by  the  Hamiltonian  function 


Notice  how  one  can  go  from  the  variational  system  (10)  to  its  adjoint  by  the 
completely  canonical  mapping 

u  =  U*,  w  =  W*,  U  =  -  u  ,  W  =  -  w  .  (19) 


This  is  a  nontrivial  example  of  the  law  of  exchange  between  coordinates  and 
momenta  in  a  phase  space. 


As  it  is  well  known,  any  orbital  integral  gives  rise  by  means  of  its 
gradient  to  a  solution  of  the  adjoint  variational  equations.  Thus,  in  view 
of  the  proposition  we  proved  in  the  preceding  paragraph,  the  gradients  VG, 
VP,  VQ  constitute  three  linearly  independent  solutions  of  the  equations  (17) 
Now  we  refer  ourselves  to  the  dual  of  Jacobi's  theorem  of  the  last  multiplier 
the  adjoint  equations  (17)  having  a  multiplier  (which  is  equal  to  1),  once  we 
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know  three  independent  solutions,  we  must  be  able  to  produce  the  general 
solution  by  means  of  quadratures  only. 


To  this  effect  let  us  introduce  four  unknown  functions  ot,8,y,6 
time  under  the  provision  that  the  linear  combinations 

* 

u  =  aC  +  6P  *-  yQ  , 

X  X  X 

w  =  aG  +  8P  +  yQ  * 

y  y  y 

U  =  aGx  +  8Px  +  yQx  +  <$» 

W*  =  aGy  +  8Py  +  yQy 


of  the 


(20) 


make  a  solution  of  the  adjoint  variational  equations  (17).  On  substituting 
the  expressions  (20)  into  the  equations  (l7),  we  find  that  the  unknowns 
a,8,Y  and  5  ought  to  be  solutions  of  the  differential  system 


dG  +  8P  +  yQ 

X  X  X 

■  -  ^(^2  -  iY- 

(21a) 

r  '  r  ' 

«Gy  +  6Py  +  YQy 

=  _  3iL2aL  6 
■^3  2  » 

r  r 

(21b) 

dGx  +  6PX  +  yQx+6 

*  o, 

(21c) 

dGy  +  8Py  +  yQy 

*  0. 

(21d) 

In  order  to  solve  it,  we  begin  by  using  the  equations  (21a)  and  (21b)  to  ex- 
press  2  and  y  by  means  of  a  and  6,  so  that 


12" 


JL 

3 

r 


3^  Q 

r 


-  - l)  % 


6  -  C12a, 


•  u 

A12Y  =  3 

r 


3^P  + 

2  x 
r 


(22a) 


y 


6  +  B12a; 


(22b) 


-li¬ 


ve  bring  these  expressions  into  the  equation  (21d)  to  express  a  in  terms 
of  6,  which  gives  the  relation 


so  that  (22c)  becomes 


a  = 


x6 

3 


Gr 


_d_ 

dt 


(¥)- 


3p 


(22c) 


(23a) 


We  then  substitute  (23a)  into  (22a)  and  (22b)  to  obtain  the  relations 


a  =  2vx  6 

B  2  3’ 
G  r 


(23b) 


.  xo 

Y'?7 


^ +  c  £(t) 


(23c) 


On  replacing  a,  B  and  y  in  (21c)  by  their  expressions  as  given  by  the 
relations  (23)  ,  we  transform  (21c)  into  the  differential  equation 
f  /  2 


a  n  1  M 

A12D1246  ~  "  3 
r 


~  1)  (A23D124"A24D123>  +  (A14D123-A13D124) 


But  from  the  definitions  (13)  and  (15),  we  compute  that 


uG 

A23D124~A24D123  *  ^3  X* 


/. 

A14D12j  ri13D124  "  'kT  Y’ 
r 
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Accordingly  che  linear  differential  equations  which  define  5  take  the 


following  simple  form 


1  _  M  _sL/_2i 


M.  _£L/.2L\ 
r3  dt(r3) 


it  is  immediately  integrable  by  quadrature  to  give  that 


5  =  K,G  — 
1  x 


(24a) 


K,  being  an  arbitrary  constant  of  integration.  When  we  substitute  this  ex¬ 
pression  of  6  into  the  equations  (23)  and  perform  the  quadratures  as  indi¬ 


cated,  we  obtain  that 


’  K2  -  K1  (3ut  - 
6  *  K3  +KlT/*dt- 

■  k4  +  K1  f  hic  +  K1  t  ■ 


(24b) 


(24c) 


(24d) 


This  concludes  the  cask  we  assigned  to  ourselves  in  this  paragraph:  in  con¬ 
formity  with  the  dual  of  Jacobi's  theorem  of  the  last  multiplier,  solving 
the  adjoint  variational  equations  (17)  reduces  to  performing  the  quadratures 
indicated  in  (24c)  and  (24d). 


3 .  The  quadratures  required  to  solve  the  adjoint  variational  equations. 

By  going  back  to  the  first  two  identities  (8) ,  we  immediately  derive 

/x  dt  =  — [G(GPt  -  ?sQr2)  -  uP  fr  dtl  , 

P  -K)  L  J  J 

fy  dt  =  —— -lorcnt  +  *i?r2)  -  uQ  /r  dtl, 


(25) 
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2  2 

provided  of  course  that  P  +  Q  f  0,  which  is  the  case  if  and  only  if 
the  Keplerian  motion  (9)  is  not  circular.  This  will  be  our  assumption  in 
the  present  paragraph.  In  view  of  (25),  the  problem  is  now  reduced  to 
constructing  a  primitive  of  r. 


Let  us  begin  by  observing  that 


(26) 


(27) 


But  Chamyi  (1963)  has  shown  how  the  two  primitives  contained  in  the  right 
member  of  (27)  are  to  be  constructed. 

For  convenience  we  put 

jih  <28) 

r  r 

On  dividing  the  orbital  integral  of  energy  written  in  the  form 

H  =  5s(^'+4)"r  (29) 

.2  .2 

first  by  r  and  then  by  r  r,  we  find  that 

2HIq  +  2ul1  -  G2I2  =  t,  (30a) 

2HI1  +  2ul2  -  G2I3  =  f&. 


(30b) 
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In  the  same  way,  on  dividing  the  equation  (26)  first  by  r  and  then  by 

,2 

r  /r,  we  obtain  that 


•‘l  -  C\  -  -  /ft  dt’ 


(30c) 


jI2  “  G  *3  / . 2  dt* 


(30d) 


However,  on  account  of  (26),  an  integration  by  parts  yields  immediately  that 


/r2dt  -  rr  -  -  *)dt. 


whence,  by  recourse  to  the  integral  (29),  we  establish  that 


ff  -  =  7  (rir  '  2Ht); 


(31a) 


this  exhibits  the  right  member  of  (30b)  in  a  closed  form.  On  the  other  hand 
we  observe  that 


iL/l\  ml.H 

dt\rj  .2 


wherefrom,  by  quadrature,  we  arrive  at  the  primitive 


Li.  .  i. 

—  dt  -  t  -  -r  ; 


(31b) 


this  produces  the  right  member  of  (30c)  in  a  closed  form.  Eventually,  because 


we  find  that 


(31c) 


and  so  we  exhibit  the  right  member  of  (30d)  in  a  closed  form.  To  conclude: 
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in  order  to  construct  t.'ie  primitives  Ig  and  1^  requested  by  (27)  ,  all 

we  have  to  do  is  to  solve  the  linear  system  (30)  in  the  unknowns  I^(n  =  0,1, 2, 3), 

provided  we  use  the  primitives  (31). 


We  use  the  equations  (30c)  and  (30 d )  to  express  1^  and  as  functions 

of  I2;  then  we  substitute  the  results  into  (30b)  and  we  find  that 

I  =  (32a) 

(P^KT)rr 

Thereafter  we  go  back  to  the  equation  (30c)  to  obtain  that 


(32b) 


(32c) 


At  this  point,  in  order  to  divide  both  sides  of  (32c)  by  H,  we  have  to 
exclude  the  parabolic  motions,  which  we  shall  do  for  the  remainder  of  this 
paragraph  and  the  next  one. 


Let  us  replace  Iq  and  1^  in  (27)  by  their  expressions  as  given  by 
(32b)  and  (32c).  After  some  cumbersome  but  elementary  transformations,  we 


finally  find  that 


/r  dt "  s  <g2  +  ur)  - 1»(™ + 


It  follows  therefrom  by  the  relations  (25)  tha’ 
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r  ,  3p 

Jx  dt 


Pr(ur  +  G  )  +  2HG0r 
4ll(P2  +  Q2) 


(34a) 


)y  dt  =  4^  ‘ 


Qr(ur  +  G  )  -  2HGPr 
4H(P2  +  Q2) 


(34b) 


After  a  few  more  manipulations  we  arrive  at  a  simpler  form  for  these  primitives 


I'  3P  12 

h  dt  =  4HC  +  4i(Xr  "  Gy)* 

/y  dc  =  4HC  +  ^(Yr2  +  Gx). 


(35a) 


(35b) 


It  should  be  said  that  we  obtain  the  right  side  of  (35a)  by  omitting  the 
3  2  2 

constant  OG  /4H(P  +Q  )  from  the  right  side  of  (34a),  and  the  right  side 

3  2  2 

of  (35b)  by  likewise  omitting  the  constant  -PG  /4H(P"+Q  )  from  the  right 
side  of  (34b). 

To  our  knowledge,  the  formulas  (33),  (35a)  and  (35b)  appear  to  be  new 
in  the  literature  of  the  problem  of  two  bodies. 


4.  General  solution  of  the  variational  equations. 


Before  we  compute  the  general  solution  of  the  variational  equations  for 
an  elliptic  or  an  hyperbolic  motion,  let  us  take  the  partial  derivatives  of 
the  basic  identity  (7)  so  that 


PP  +  QQ  =  G(uG  ~  +  2HY), 
x  x  3 


PP  +  QQ  =  G(uG  -*r  -  2HX)  , 

y  y  r3 

PPV  +  QQ..  =  G(  Gx  -  2Hy)  , 

X  A 


PPy  +  QQy  =  C(  Gy  +  2Hx) . 
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In  view  of  these  formulas  we  draw  from  the  quadratures  (35)  that 
pxix  dt  +  Qx/y  dt  =  m 

pyix  dt  +  %fy  dt  =  m 

pxix  dc  +  ojy dc  =  4 

py/x  dt  +  qjydt  =~ 

We  are  now  ready  to  substitute  the  functions  a,  $,  y,  6  given  by  the 
formulae  (24)  into  the  linear  combinations  (20).  A  few  trivial  manipulations 
yield  the  final  resulL: 

u*  =  K'(X  -  Sir^t)  +  K2Gx  +  K3Px  +  KQ 
r 

v*  =  K!  (Y  -  3u*4t)  +  K„G  +  K.P  +  K.Q  , 

1  J  2  y  3  y  4  y 

r  (38) 


3t 

+ 

2HY] 

-  2Hy  - 

GX 

'  r 

3t 

(u° 

- 

2HX 

+  2Hx  - 

GY 

. 

V  r 

7 

V 

3t 

(  GX 

- 

2Hyy 

-  2G2x 

* 

3t 

^GY 

+ 

2Hx 

1  -  2G2y 

* 

(37) 


U*  =  K’(2x  -  3Xt)  +  K2Gx  +  K3Px  +  K4Qx, 

V*  =  K^(2y  -  3Yt)  +  +  K^, 

it  being  understood  that  the  arbitrary  constant  of  integration  is  now 

replaced  by  a  constant  such  that  2HK3  =  -yGK^. 

We  find  it  convenient  to  define  the  vector 

S  =  s1&1  +  +  s3&3  +  s4&4 


whose  components  are 
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S1  =  X  -  3u  ~  t, 
r 


S3  =  2x  -  3Xt, 


S.  =  Y  -  3p  t,  S.  =  2y  -  3Yt. 


We  just  proved  that  S  is  a  solution  of  the  adjoint  variational  equations 
(17). 

Let  us  compute  the  one-dimensional  4-vector  £  a  VG  a  VP  a  VQ;  we  find 


S  a  A  7P  A  7Q  =  -2HG  e,  a  a  e->  a  e,  . 


This  proves  that  the  matrix 


* 

A  (c) 


whose  columns  are  the  vectors  jS,  VG,  VP,  VQ  taken  in  that  order  has  its 

2 

determinant  equal  to  -2HG  ;  it  thus  constitutes  a  fundamental  set  of  solu¬ 
tions  to  the  adjoint  variational  equations. 


In  order  to  compute  the  inverse  B  (t)  of  A  (t) ,  we  compute  the  fol¬ 
lowing  3-vectors: 


s  AVG  A  VP  =  Z  E  £.A  &  A  v 
(i*j,k)  J 

S  a  VG  a  VQ  =  A  £j  A 

S  A  VP  A  VQ  =  ^  k)  Jijk  ^i^jA  V 


the  components  are 


E 1 2 3  -  (Yr  +  Gx  +  3Qt) 

r 
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E124  =  (Yr2  +  Gx  +  3Qt)  -  2HG, 

r 

E134  =  -  Y  (Yr2  +  Gx  +  3Qr)  +  2Hy2, 

e234  =  x  ^Yr2  +  Gx  +  3Qc>  ~  2Hxy» 

F12 -i  =  (Xr2  -  Gy  +  3Pt)  +  2HG, 

r 

Fi24 =  nr (Xr2  ■ Gy  +  3Pt)  * 

r 

Fl34  -  Y  (Xr2  -  Gy  +  3Pt)  -  2Hxy, 

F234  =  -  X  (Xr2  -  Gy  +  3Pt)  +  2Hx2, 

J123  =  2HG  (Y  -  3ir^t)  , 
r 

Ji24 =  -2HG  (x  ■  3nt) » 
r 

J134  =  -2KG  (2y  -  3Yt) , 

J234  =  2HG  (2x  -  3Xt). 


But,  as  we  know  from  Linear  Algebra,  these  3-vectors  and  the  one  we  computed 

*  * 

in  (17)  yield  the  elements  of  the  inverse  matrix  B  from  the  following 

relations : 


2 

* 

2 

* 

-2hg: 

bll  =  °234  ’ 

-2HG7 

b21 

1 

u, 

to 

OJ 

*• 

-2HG2 

b12  “  -D134 ’ 

-2HG2 

* 

boo  = 

W*. 

J134» 

2 

* 

„  2 

* 

-2HG7 

b13  =  °124 ’ 

-2HG7 

b23 

_J124’ 

2 

•k 

O 

* 

-2HG7 

b14  =  ~U123’ 

-2HG7 

b24  " 

J123* 

0 
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-2HG? 

* 

b31  = 

F 

234’ 

-2HG^ 

* 

b41  = 

-E234  * 

2 

ft 

„  2 

* 

-2HG7 

b32  “ 

-F134  * 

-2HG7 

b42  = 

E134* 

2 

* 

2 

* 

-2HG7 

b33  = 

FlzV 

-2HG7 

b43  * 

~E124* 

2 

■k 

2 

* 

-2HG7 

b34  = 

-F 

r123> 

-2HG7 

i! 

E123  ’ 

The  matrix  A  (t)  being  a  fundamental  set  of  solutions  to  the  equations 
(17) ,  the  matrix 

R*(t;t0)  =  A*(t)o  B*(t0)  (39) 

is  also  such  a  set.  But  obviously 

R*(to;to)  =  A*(co)o  B*(to)  =  V 

Hence  R  ( t: ; t^)  is  the  resolvent  of  the  adjoint  variational  equations. 


is  the  resolvent  of  the  variational  equations.  A  simple  expression  of 


R(t;tg)  is  obtained  by  putting 


A(c)  =  J»A  (t). 


(41a) 


B(t)  *  B*(t) °J_1 


(41b) 


so  chat 


R(t;t0)  =  A(t)«>B(t0). 


(42) 


The  elements  of  the  matrices  A(t)  and  B(t)  are  presented  in  Tables  1  and 


Table  I  -  The  factor  A(t)  of  the  matrizant  R(t,tg) 


all = 

2x 

-  3Xt 

a12  = 

-y 

a21  = 

2y 

-  3Yt 

a22  = 

X 

a31  = 

-  X 

+  3vr^t 
r 

a32  = 

-i 

a41  = 

-  Y 

+  3^t 

y 

a42  = 

X 

13  = 

al4 

23  =  + 

2xY 

a24 

2 

=  -Y^  + 
33 

a34 

r 

43  XY' 

a44 

-xY  +  2yX 
-xX 


XY-  ^ 
r 

2  X2 
-X  +  u~ 
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Table  II  -  The  factor  B(t)  of  the  matrizant  R(t;tg) 


5 .  The  matrizant  in  the  orbital  frame  of  reference. 

The  resolvent  (42)  is  related  to  inertial  Cartesian  frames  of  reference 

in  the  original  space:  given  a  displacement  ^-(AXq,  AyQ,  AxQ,  AyQ)  of  the 

initial  position  (xQ,y0)  and  of  the  initial  velocity  (x0,yQ)  with  respect 

to  that  fixed  frame,  it  provides  in  the  linear  approximation  the  evolution 

A(t)  * (Ax(t) ,  Ay(t) ,  Ax(t) ,  Ay(t))  of  that  displacement  along  the  orbit  T, 
/■>* 

namely 

A(t)  «  R(t;t0).A0. 
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But  there  are  other  frames  of  reference,  not  fixed  in  space,  to  which 
it  is  sometimes  more  convenient  in  practice  to  refer  the  Keplerian  orbit. 

In  this  paragraph,  we  shall  consider  the  orbital  frame  of  reference:  the  origin 
0  is  the  center  of  attraction,  the  axis  Ox*  is  the  line  from  0  to  the  posi¬ 
tion  occupied  by  the  particle,  and  the  axis  Oy'  is  taken  orthogonally  to  the 
axis  Ox'  and  90°  forward  of  the  axis  Ox'  in  the  sense  of  the  motion  of 
the  particle  on  its  orbit. 


If  6  is  the  azimuth  of  the  particle  in  the  inertial  frame  of  reference. 


such  that 


x  =  r  cos  0,  y  **  r  sin  0, 


(43) 


then  its  Cartesian  coordinates  in  the  orbital  frame  Ox'y'  are 


x'  =  r,  y'  *  0, 

and  the  Cartesian  components  in  the  frame  Ox'y'  of  its  velocity  with  respect 
to  the  frame  Oxy  are 

X'  *  r,  Y»  -  r0  -  £  .  (44) 

Let  (u1,  w* ,  U' ,  W')  be  the  components  in  the  frame  Ox'y*  of  the  displace¬ 
ment  u  in  the  frame  Oxy;  obviously 


u 

-  u' 

cos 

0  +  w' 

sin 

w 

*  — u ' 

sin 

0  +  w' 

cos 

e. 

U 

=  U' 

cos 

0  +  W 

sin 

0, 

W 

=  -U’ 

sin 

0  +  W 

cos 

0. 

(45) 


Thus  the  components  u' ,  w' ,  U' ,  W'  of  the  displacement  can  be  interpreted 


as  follows: 


0 


-24- 


(i)  u'  is  Che  correction  in  radial  position  or  range; 

(ii)  w*  the  correction  in  transversal  position; 

(iii)  U'  the  correction  in  radial  velocity  or  range  rate; 

(iv)  W*  the  correction  in  transversal  velocity. 

Notice  that  the  relations  (45)  define  a  time-dependent  canonical  tranformation ; 
it  can  be  taken  as  being  generated  by  the  function 

S'  =  S'(U,W,u',w’)  =  -U(u'  cos  0  +  w'  sin  0)  -  W(-u'  sin  0  +  w'  cos  0) 

so  that  its  remainder  is  the  function 

!£'«  -  (uW  -  wU)0  «  -G(uW  -  wU) 0.  (46) 

O  t 

The  transformed  variational  Hamiltonian  is  the  function 

V  =  V  (u'.w’.U'.W’)  «<v  +  |£'.  (47) 

d  t 

Since  the  transformation  (45)  is  a  proper  rotation  in  the  plane  of  configu¬ 
rations  and  in  the  plane  of  moments,  we  have  immediately  that 

U2  +  W2  =  U’2  +  W’2,  (48a) 

uW  -  wU  =  u’W  -  w’U’;  (48b) 

we  also  compute  that 

2  /  2 

“  1^  u2  +  6xy  +  ^3^2  -  ljw2  =  2u’2  -  w’2.  (48c) 

This  last  expression  is  interesting  in  that  it  exhibits  that  the  orbital  axes 
Ox’  and  Oy’  are  the  principal  axes  of  that  quadratic  form  which  is  the 
force  function  in  the  variational  Hamiltonian  'V.  Thus,  mathematically 
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speaking,  the  orbital  axes  Ox*  and  Oy'  are  introduced  to  diagonalize 
partially  the  variational  Hamiltonian. 


On  introducing  (48a),  (48b)  and  (48c)  into  (11)  and  (46),  we  eventually 
obtain  from  (47)  that 

V  =  5i(U,2+W'2)  -  ~  (u'W-w’U1)  -  ~  (2u,2-w’2)  (49) 

r  2r 


The  variational  equations  derived  from  'V*  are 


=  _  u'  +  —  w'  U»  =  _  —  =  JL  y'  +  JL  u» 

u  ,  u  -r  9  w  ,  u  ,  w  -r  u  , 


3u*  2 

r 


3V  G  ,  •’  W  G  u  , 

”  '  3W '  ‘  "  -  ~  “  '  "  ’  -  St’  *  -  “2  U  -  7  j" 

r  r  2r 


in  their  Hamiltonian  form,  or 


u*  -  2~  w' 
r 


G  +ur  ,  .  „Gr  , 

- 7—  u  +  2~  w' 

4  3 

r  r 


w'  +  2~  u* 


Gr  ,  G  -ur  , 
3  u  -  4  w 

r  r 


in  their  Lagrangian  form.  These  are  the  equations  which  Brumberg  (1961)  and 
Charnyi  (1964)  considered  in  preference  to  the  equations  (10) . 


We  are  able  to  construct  in  a  direct  way  the  matrizant  R’(t;tg)  of 
the  system  (50) .  Let  us  put 


I  cos  0  sin  0 
-sin  6  cos  0 


cos  0  sin  0 

-sin  0  cos  0 


.  L’ 


cos  0  -sin  0 
sin  0  cos  0 


cos  0  -sin  0 
sin  0  cos  0, 


-26- 


because 


\ 

/  f  \ 

/  f  \ 

1 

u  1 

i  i 

f  U 

f  u  1 

1  1 

1  u 

w 

11 

w’ 

> 

w' 

II 

w 

U 

|U,| 

U’ 

|  | 

U 

I 

\  Wi 

\w  i 

IwJ 

\  W 

we  have  that 

R’(t;t0)  -  K' oR(t ;tQ) °L' . 

This  suggests  introducing  the  matrices 

A'(t)  =  K'oA(t) ,  B'(t)  -  B(t)°L',  (52) 

so  that  the  matrizant  is  decomposed  into  the  product 

R'(t;t0)  =  A'COoB'Uq).  (53) 

Tables  III  and  IV  present  the  elements  of  the  factor  matrices  A'(t)  and 
B ' ( t)  which  compose  the  matrizant  of  the  elliptic  and  hyperbolic  Keplerian 
motions  (G  ^  0,  e  j4  0,1)  in  the  orbital  frame  of  reference. 

Table  III  -  The  factor  A'(t)  of  the  matrizant  R* (tjt^) 
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6.  Hill 1 s  equation  for  the  problem  of  two  bodies. 

Because  we  restrict  ourselves  to  elliptic  and  hyperbolic  motions, 
nowhere  along  the  orbit  does  the  velocity  vanish.  Hence  we  define  without 
ambiguity  an  angle  4>  such  that 

X  *  V  cos  Y  »  V  sin  <fi;  (54) 

4>  is  the  angle  which  the  velocity  vector  makes  with  the  inertial  axis  Ox. 

From  the  definition  (54) ,  we  immediately  draw  that 

XX  +  YY  -  VV,  XX  +  YY  «  VV  -  V2*2, 

.  2'  ....  ..  2.. 

XY  -  XY  *  VI,  XY  -  YX  -  2W$  +  V  4>. 

Now  let  us  for  a  moment  consider  a  more  general  dynamic  system  described 
by  the  Hamiltonian  function 

* s(X2  +  Y2)  -OKx.y).  (56) 

On  using  the  differential  equations  of  motion  derived  from  (56)  as  well  as 
the  identities  (55)  ,  we  obtain  that 

or  cos  4>  +  0T  sin  $  «  V, 
a  y 

(57) 

OJr  sin  -OT  cos  <f>  -  -V$, 

a  y 

^xx  COs2  ♦  +  2f>ky  COS  ♦  Sln  ♦  +  f*yy  sin2<(’  “  V  ■  ^2' 

Orxy(cos  $  -  sin  $)  +  <0 ryy  -  or^)  cos  sin  <j>  -  2-^  4>  +  $, 

°xx  Si°2  *  "  C0S  ♦  sin  ♦  +0^y  cos2*  m%x  +  (*yy  +  ^  ~  V* 
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We  transform  the  displacement  (u,w,U,W)  in  the  following  way: 


u  ■  u"  cos  4>  -  w"  sin  4> , 
w  =  u"  sin  d>  +  w"  cos  <p, 


U  =  U"  cos  <p  -  W"  sin  4> , 
W  =  U"  sin  <p  +  W"  cos  <p. 


(58) 


The  components  u",  w",  U",  W"  of  the  displacement  can  be  interpreted  as 
follows : 

(i)  u"  is  the  tangential  variation  or  the  variation  along 
the  track; 

(ii)  w",  the  normal  variation,  or  the  variation  across  the  track; 

(iii)  U",  the  correction  in  tangential  velocity; 

(iv)  W",  the  correction  in  normal  velocity. 


The  relations  (58)  define  a  time-dependent  canonical  transformation;  it  can 
be  considered  as  being  generated  by  the  function 

S"  -  s"(U,W,u",w")  *  -U(u"cos  +  w"sin  <J>)-W(-u"sin  $  +  w"eos  4>) 


so  that  its  remainder  is  the  function 

*  -  (uW  -  wU) 

But  the  variational  Hamiltonian  associated  with  (56)  is  the  function 

'V  **  h(U2  +  W2)  -  h(0^x u2  f  20£yiiv  +Ofyyv2) 

which,  in  virtue  of  the  canonical  transformation  (58),  transforms  into  the 
function 

V  5  V,(u",w",U",W")  -  -V+  |f  . 

d  Z 

In  view  of  the  transformation  relations  (58)  and  of  r^_  identities  (57), 


we  readily  obtain  that 
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V  -  is(U"2  +  W,,2)-$(u"W"  -  w"U") 


(59) 


-h 


[(i  -  i2)u"2  +  2(4  *  +  i)»v  +(orxit  +  o$,y  +  *2  - |)»''2] 


From  the  Hamiltonian  equations 


u" 


u"  +  <j>w", 

w"  -  W"  -  $u", 

U"  -  4>W"  +  (|  -  u"  +  {jjfy  +  $)w", 

w"  -  ~4>u"  +  (2^  +  $)u"  +  (orxx  +o^y  +  V 


-  v)*"’ 


we  deduce  the  Lagrangian  equation  for  the  normal  variations 

w"  +  -  2$2  -»xx  “Ofyyjw"  -  2  |  (Vu"  -  Vu"). 

But  the  variational  integral  of  energy 


(60) 


(61) 


<$H  -  XU  +  YV  -  Ofxu  -  OJryv 


becomes  in  the  intrinsic  variations  the  function 


6H  -  VU"  -  Vu"  -  V$w" 


(62) 


or,  on  using  the  second  of  the  equations  (60), 


<5H  -  (Vu"  -  Vu")  -  2V<J>w" 


(63) 


Therefore  the  Lagrangian  equation  (61)  transforms  into  the  second  order 
linear  differential  equation: 


«"  +  0w"  -  -2^5H 


(64) 
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whose  coefficient  is  the  function 


V  •  2 

9=-+24>  -or  -or  . 

V  xx  yy 


(65) 


To  conclude,  the  resolution  of  the  Hamiltonian  system  (60)  reduces  to  the 
integration  of  a  second  order  differential  equation  in  the  normal  displace¬ 
ment  w  and  the  quadrature 


—  -5i  a  +  2-^w" 

dt  v  \r 


(66) 


as  indicated  by  the  variational  integral  (63)  for  the  tangential  displacement. 

We  now  return  to  the  problem  of  two  bodies.  In  this  case,  in  view  of 
the  identities  (55)  and  of  the  orbital  equations  of  motion,  we  find  that 


V  =  - 


J iL  « 

v  2 

Vr 


UG 

u2  3 

V  r 


V  -  -  r)] 

Vr  L  V  r  V  r/J 


°*xx  +%y  *  ^ 


Accordingly,  the  coefficient  (65)  in  Hill’s  equation  (64)  is  the  function 

,2 


e 


V2r 


h  [»*  -  ^ 2  -  *)] 


(67) 


which  can  also  be  written  as 


0 


„2  3 

V  r 


jv2  +  3G($>  -  &)J 


to  provide  a  better  indication  of  its  physical  meaning. 
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7.  The  matrizant  In  the  intrinsic  frame  of  reference. 


In  the  same  way  as  we  did  in  the  orbital  frame  of  reference,  we  decom¬ 
pose  the  matrizant  R"(t;tQ)  in  the  intrinsic  frame  of  reference  into  the 
product 

R"(t;t0)  *  A"(t)oB"(t0), 

with 

A"(t)  =  K"°A(t) ,  B"(t)  -  B(t) °L", 


1  cos 

4> 

sin  $ 

0 

0 

1 

1  cos 

4> 

sin  4> 

0 

0 

1 

K”  = 

-sin 

cos  4> 

0 

0 

9 

L"  = 

-sin 

♦ 

cos 

0 

0 

0 

i 

0 

cos  <P 

sin 

*1 

0 

0 

cos  $ 

sin 

* 

\  0 

0 

-sin  <p 

cos 

*1 

1  0 

0 

-sin  $ 

cos 

♦/ 

The  elements  of  the  matrices  A"(t)  and  B"U)  are  listed  in  Tables  V  and  VI. 

We  denote  by  a|'  the  vector  whose  components  are  the  elements  of  the 
tbi 

i  column  in  the  matrix  A(t).  For  1  <  i  <  4,  the  vector  a",  is  a  solu- 

—  **  ~  i 

tion  of  the  variational  equations.  Let  us  compute  the  value  taken  by  the 
variational  integral  of  energy  along  that  variation;  we  find  that 


(^la’V  *  -2H»  and  (VH|a'y  -  0  (2  £  j  £  4). 


(68) 


Variations  such  that  (VH|u)  ■  0  are  called  isoenergetic.  Thus  what  the 
equalities  (68)  prove  is  that  the  isoenergetic  variations  constitute  a  three- 
dimensional  space,  and  that  the  variations  a"^  (j  »  2,3,4)  form  a  base  of 


that  vector  subspace. 
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Table  VI  -  The  factor  B"(t)  in  the  matrizant  R"(t;tQ) 


h"  -  -Ji¬ 
ll  “  2H 


b» 

12  2H 


h _  V  JT 
b21  ~  ~  G  '  ^ 


b22  '  -  ❖ 


b31  "  ^2(Xr2  ‘  <*  +  3Pt)  +  S 


b32  -  -^(Xr2  -  Gy  +  3Pt)  +  ~ 
2HG 


b41  *  ~^~2^r2  +  Gx  +  3Qt)  -  ~ 
2HG 


b"  =  ~^r(Yr2  +  Gx  +  3Qt)  +  ^ 
2HG 


h"  -  -Si- 
13  2H 


b"  »  0 
14 


b..  »  2rr  _  ^ 

23  GV 


b"  =  -  - 

24  V 


b"  -  — ^-r(Xr2  -  Gy  +  3Pt)  +  2~ 
2HG  VGZ 


h"  «  JL 
34  "  GV 


b"  *  -  — ^(Yr2  +  Gx  +  3Qt)  + 

2HG  VG 


b"  -  -  -St- 
44  GV 
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Fore  most  among  the  isoenergetic  variations  is  the  solution 

&o  -  <vvW  such  that 


un  -  X,  wn  =  Y,  Un  -  X,  Vn  «  Y; 


its  components  in  the  intrinsic  frame  of  reference  are 


u"  -  V,  v"  =  0,  U"  -  V,  W"  =  -V* 


The  reader  could  check  that  this  remarkable  isoenergetic  displacement  is  the 
linear  combination 


•^0  ~  +  60&3  +  Y0% 


whose  coefficients  are  the  functions 


a0  "  “  4*  “  2t* 


s0  *  i(v2  -  7) +  ^  + • 

U 

Y0  *  i(y2  -  r)  +  $  +  TCI  • 


One  usually  completes  the  base  of  the  vector  subspace  of  isoenergetic  vari¬ 
ations  by  adjoining  to  the  vector  (70)  the  solutions  ^  and  u^  defined 
by  the  following  initial  conditions: 


u'^vtg)  »  0,  wi<to)  “  x*  w^(t0)  =  0, 

u^(tQ)  «  0,  w^(tQ)  "  °»  ,J2(t0)  *  1* 


The  variation  results  from  displacing  the  initial  position  along  the 

normal  to  the  orbit  at  that  point  without  modifying  the  velocity,  while  the 


0 
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variation  ^  is  caused  by  an  initial  change  of  the  normal  component  of 
the  velocity  (or  for  that  matter  an  initial  impulse  in  the  normal  direction) 
without  changing  the  initial  position.  To  make  such  displacements  isoenergetic, 
we  adjust  the  change  in  the  tangential  component  of  the  velocity  according  to 
the  variational  integral  of  energy  (63)  in  which  we  put  6H  -  0.  There  re¬ 
sults  that 

uJ(tQ)  -  2*(t0),  uJ(tQ)  =  0.  (72) 

Since  we  assume  that  they  are  isoenergetic,  the  variations  and  ^  are 

linear  combinations  of  the  kind 

%  “  al*S  +  61*3  +  Yl*4’ 

~2  =  “2^  +  +  V*  ’ 


They  satisfy  the  initial  conditions  (71)  if  and  only  if 


Accordingly  their  components  are  as  follows 

“I  •  +  <Vo-Vo>T  -  (Vo-Vo>4 

ui  ■  -Vto*  +  (vovVo)(!cxtv5|  +  <voVVo)(2gy-v2*)]  • 

GVov 

"i  ■  kvVo'«  -  <v<fW*]  • 

Gvo 
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W1  =  [2GV  +  (V0X0~Vo)(VY+^y)  “  (V0V0'Vo)(VX^x)]  (73) 

Gvo L 

for  the  basic  isoenergetic  displacement  and 

u2  =  V^(G  -  YX0  +  V  * 

W2  =  G^?[2Gri'  ‘  <v2y+2GX>xo  +  (v2x"2GY)y01, 

.°  (74) 

U2  “CT0(yX0  ‘  V’ 

w2  *  gv^[2GV  ■  (VY+^y)xo  +  (VX+Vx)y0] 

for  the  other  basic  isoenergetic  displacement  u^. 
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Conclusion 

The  variational  equations  of  the  problem  of  two  bodies  have  been  solved 
from  first  principles  applying  Jacobi's  dual  theorem  of  the  last  multiplier 
to  the  adjoint  of  the  Hamiltonian  system  expressed  in  Cartesian  coordinates 
with  respect  to  an  arbitrarily  given  inertial  frame  of  reference.  Once  the 
matrizant  R(t;tQ)  in  that  system  has  been  obtained  in  closed  form,  other 
forms  of  the  resolvent  are  derived  simply  by  homogeneous  canonical  extensions. 
By  way  of  illustration,  we  have  recovered  in  this  way  the  matrizants  in  the 
orbital  frame  of  reference  and  in  Hill's  intrinsic  coordinate  system.  To 
our  knowledge  this  last  resolvent  is  not  to  be  found  elsewhere  in  the  litera¬ 
ture  on  the  problem  of  two  bodies. 

The  treatment  we  give  in  this  paper  is  incomplete  in  many  respects.  For 
instance  some  who  deal  with  orbital  transfers  or  trajectory  designs  by  matching 
conics  would  like  to  introduce  a  so-called  universal  variable  which,  for 
Keplerian  motions  as  functions  of  the  eccentricity  e,  acts  as  a  locally 
uniformizing  variable  (Sconzo  1967)  in  the  neighborhood  of  e  *  1.  Goodyear 
(1965)  has  been  the  first  to  introduce  a  universal  time  variable  in  the  matri¬ 
zant  of  the  problfc-i  of  two  bodies;  he  did  it  as  an  extension  of  the  method 
proposed  by  Stumpff  (1947,  1959,  1962)  for  the  orbital  equations.  At  this 
point,  the  problem  is  to  find  a  w~y  of  simplifying  the  lengthy  and,  at  places, 
cumbersome  manipulations  performed  by  Goodyear  (1966). 
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